Orbital-free energy functional for electrons in two dimensions 
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We derive a non-empirical, orbital-free density functional for the total energy of interacting elec- 
trons in two dimensions. The functional consists of a local formula for the interaction energy, where 
we follow the lines introduced by Parr for three-dimensional systems [R. G. Parr, J. Phys. Chem. 
92, 3060 (1988)], and the Thomas-Fermi approximation for the kinetic energy. The freedom from 
orbitals and from the Hartree integral makes the proposed approximation numerically highly effi- 
cient. The total energies obtained for confined two-dimensional systems are in a good agreement 
with the standard local-density approximation within density-functional theory, and considerably 
more accurate than the Thomas-Fermi approximation. 
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I. INTRODUCTION 

Two-dimensional (2D) electronic systems have at- 
tracted vast interest since the beginning of semiconduc- 
tor technology. Important examples are quantum-Hall 
systems and different types of quantum-dot (QD) de- 
vices i Technological development has also increased the 
need for computational methods capable to deal with the 
many-electron problem in reduced dimensions. Among 
the available methods is the well-known local-density 
approximation (LDA) within the celebrated density- 
functional theory^ (DFT). The 2D-LDA consists of the 
exchange functional derived for the homogeneous 2D 
electron gas^ and the corresponding correlation func- 
tional constructed using quantum Monte Carlo meth- 
ods^ At present, DFT with the 2D-LDA, and es- 
pecially their spin-dependent (and current-dependent) 
extensions, are among the standard methods in the 
electronic-structure calculations of semiconductor QD's. 6 
Further developments of 2D density functionals have be- 
gun very recently for both exchange 7 and correlation. 8 ' 9 

Although the LDA, for example, is an explicit den- 
sity functional, so that the total density is the sole input 
variable (instead of the electronic orbitals), the standard 
Kohn-Sham (KS) scheme in DFT still requires the com- 
putation of the KS orbitals for the single-particle kinetic 
energy. This sets limitations to the number of electrons 
that can be treated numerically. The so-called orbital- 
free DF T 10 ' 11 : 12 scales better in this respect, but within 
this approach the construction of an accurate energy 
functionals (in particular for 3D systems) has resulted 
to be a complicated task. The "traditional" Thomas- 
Fermi (TF) approximation may serve as an important 
example of an orbital-free functional. The TF approach 
has been put on a mathematically rigorous basisji 3 - and 
also analyzed in 2D in detail by Lieb et alJ^- Further- 
more, the TF theory has been successfully applied in 
the electronic-structure calculations of, e.g., quantum- 
Hall systems, where the importance of e-e interactions 
has been addressed^ The TF energy functional has, 



however, the obvious deficiency to treat the e-e inter- 
action only classically (i.e., only Hartree energy in in- 
cluded). Therefore, in the regime of small number of 
particles and/or low densities (strong interactions) the 
performance of the TF method is highly questionable due 
to the lack of quantum mechanical effects (exchange and 
correlation) . 

In this paper we aim at bridging the gap between the 
numerical efficiency of TF method, and the accuracy of 
standard KS-DFT, for electronic-structure calculations 
in 2D. To this end, we present an explicit density func- 
tional for the total energy which accounts for the classi- 
cal and, to some extent, also for the quantum mechanical 
contribution to the interaction energy. In the derivation 
for the interaction energy we follow the general lines al- 
ready employed in the 3D case by ParrJ£ In particular, 
we apply a Gaussian approximation of the second-order 
density matrh* 8 - and make use of the properties of the 
interaction energy under a scaling transformation. Com- 
bining the resulting formula with the TF approximation 
for the kinetic energy leads to an explicit density func- 
tional for the total energy. Applications to 2D QD's and 
rectangular quantum slabs (QS's) up to hundreds of elec- 
trons show a significant improvement over the TF ener- 
gies when compared with the LDA results in a wide range 
of the e-e interaction strength. 



II. DERIVATION OF THE APPROXIMATION 

Our aim is to obtain a good estimation of the total 
energy of a 2D system with a large number of electrons 
using a simple and computationally convenient formula. 

First, let us consider the the e-e interaction energy. 
This can be expressed in terms of the spinless second- 
order density matrix as 

W= fdrj dr 2 ^4 (1) 
J J r x -r 2 
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where 



P2(ri,r 2 ) 



N(N - 1) 



dN 



x |*(n£7i,r 2 <j 2 ,3,...,iV)| 
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Here, 2, N) stands for the ground-state many- 
body wave function and J dN denotes the spatial integra- 
tion and spin summation over the TVth spatial spin coor- 
dinate (r^vcrjv)- Hartree atomic units are used through- 
out the paper unless stated otherwise. The above defini- 
tion implies the normalization 



N(N-l) 



J dri J dr 2 P2(ri,r 2 ), 



(3) 



and pi can be interpreted as the distribution density of 
the electronic pairs. 

Next, we will specialize all the expressions to the 2D 
case and derive a local-density approximation for the in- 
teraction energy W defined in Eq. In the average, 
r = (ri + r 2 )/2, and relative, s = ri — r 2 , coordinates, 
Eq. ([T]) can be rewritten as 



W = 2ir J dr J ds p 2 {r,s), 



(4) 



where we have introduced the cylindrical average of p 2 , 
which is defined as 



P2(r,s) = ^-^ d<f> s p 2 ( 



s s 
l+ 2' r -2 



(5) 



We assume a Gaussian approximation to be valid for the 
cylindrical average of the pair-density distribution func- 
tion, 



p 2 (r, s) » p 2 (r, r)exp 



(6) 



where we have introduced /3 2 (r) as a quantity to be de- 
termined below. Substituting Eq. © in Eq. and 
integrating over the relative coordinate, we obtain 

W^ 3 / 2 |drp 2 (r,r)/3 2 1/2 (r). (7) 

Similarly, substituting Eq. ([6]) in Eq. Q we obtain 

N(N-1) = 2tt /drp 2 (r,r)/3 2 (r). (8) 



An additional and crucial assumption is introduced by 
imposing the integrands of Eqs. ([7|) and JHJ), respectively, 
to be dependent on the space variable through the par- 
ticle density. Thus, we may write 



and 



p 2 (r,r) = p 2 (p(r)), 
/? 2 (r) =&(/>«)■ 



(9) 
(10) 



It is possible to work out the dependencies on the particle 
densities of the above quantities by a dimensional argu- 
ment. Under uniform scaling of the coordinates, r — > Ar 
(with < A < oo), and of the norm-preserving many- 
body wavefunction 



*x(t u ...,t n ) = A JV *(Ar 1 ,...,Ar 



(11) 



As a consequence, the other quantities of interest scale 
as 



and 



P2,A(ri,r 2 ) = A 4 p 2 (Ari, Ar 2 ), 



p x (r) = A 2 p(Ar), 



W[$x] = XW[^}. 



(12) 



(13) 



(14) 



By using the assumptions in Eqs. ([9]) and (fT0|) together 
with the scaling properties listed above, and by a dimen- 
sional argument, we arrive at the following expressions 
for the integrands in Eqs. J7J) and ([5]), respectively: 



and 



p 2 (r,r)/3 2 1/3 (r)=C lP 3 / 2 (r), 



p 2 (r,r)/3 2 (r)=C 2 p(r), 



(15) 



(16) 



where Ci, and C 2 are constants. Eqs. (fl5|> and (TT 
imply 



c\ 



p- 1 (r)=Ap-\r), 



and 



P2 (r,r) = ^p 2 (r) = 5p 2 (r). 



(17) 



(18) 



An estimation for the latter factor B can be obtained by 
considering the Hartree-Fock (HF) case, for which 



p 2 ,HF(r,r) = ipH F (r). 



(19) 



Hence, B = 1/4. The other factor A can be determined 
by imposing the normalization condition in Eq. |8|). 



A = 



2{N - 1) 



(20) 



Now we have all the information to give an explicit ex- 
pression for the interaction energy, which results to be 



W[ P ] = y^ I drp^r). 



(21) 



We emphasize that the expression gives W = for N — 
1, while this is not recovered by the LDA and the TF 
approximations. Of course, for a large electron number 



3 



(N 3> 1) the above expression can be simplified as N — 

An interesting feature in Eq. (|2ip is the fact that it 
allows us to approximate the total e-e interaction in a 
very simple fashion, which is computationally appealing 
for systems with a large number of electrons. However, 
some caution is in order: In the derivation above, we have 
introduced the assumption in Eq. §§§ and then invoked 
the HF case in determining the coefficient B in Eq. (|T9|) . 
Alternatively, one may refer to the exact exchange ap- 
proximation within DFT. In any case, Eq. (TH)]) is valid 
for a pair density matrix coming from a wavefunction 
of the form of a single Slater determinant. As a conse- 
quence, the resulting approximation may result to be bi- 
ased toward the Hartree plus exchange energy. Neverthe- 
less, we make this choice for methodological simplicity. 
Moreover, as shown below, the resulting approximation 
allows to deal with strongly interacting systems to a very 
good extent. Alternative choice for B could be made by 
considering correlated pair densities of reference systems, 
such as the homogeneous 2D electron gas^ or by using a 
coupling-constant average wich allows to account for the 
correlation contribution to the kinetic energy^ 

We also point out that Eq. ([2Tj) has the disadvantage 
to define a functional that is not size-consistent. In fact, 
because of its nonlinear dependence on the particle num- 
ber N, even in the case that the exact kinetic energy 
would be known, the total energy of two non-interacting 
fragments is not equal to the sum of the two fragment 
energies calculated separately. 

Now, as a simple approximation for the many-body 
kinetic energy we propose the TF expression 

TtfM = | J drp 2 (r). (22) 

First of all, it is reassuring to see that the TF kinetic en- 
ergy scales as the exact one. In fact, from Eqs. (JTTJ) and 
([13]) it is straightforward to verify that T[W A ] = A 2 T[*], 
and Ttf[px] = X 2 T T f[p\- Moreover, in 2D the TF kinetic 
energy functional is particularly attractive, since its gra- 
dient corrections vanish to all orders ; 19 ' 20 ' 21 ' 22 whereas 
in 3D the first fi 2 order correction is the well-known von 
Weizsacker correction^ Besides, for the 2D Fermi gas in 
a harmonic trap the TF kinetic energy yields the exact 
non-interacting kinetic energy when the exact density is 
used as the inputs But in interacting systems, even in 
the best case the present approximation misses the corre- 
lation contribution to the kinetic energy. As mentioned 
above, it could be possible to account for this contribu- 
tion by introducing a coupling-constant average, which 
is, however, beyond the scope of this work. 

Combination of Eq. (f2"Tj) with Eq. ([2"2"|) yields an 
orbital-free density functional for the total energy, 

E[p] = T TFW + |^Ei J dvp^{v) 

+ / <Zrp(r)y ext (r). (23) 



We remind that the standard TF approximation for 
the total energy is given by 

£tf [p] = Ttf [p] + E H [p] + J dr p(r) u ex t (r) , (24) 
where the Hartree energy is defined by 

In the LDA the total energy has an expression 

E LDA [p] = T KS [p] + E H [p] + E^ A [p] + £ c LDA [p] 

+ J drp{r)v ext (r), (26) 

where the KS kinetic energy is calculated from the KS 
orbitals. Therefore, the LDA expression is an implicit 
density functional in contrast with the explicit density 
functional in Eqs. and 



III. NUMERICAL PROCEDURE 

We test our total-energy functional given in Eq. (|23[) on 
parabolic (harmonic) QD's and rectangular QS's, respec- 
tively. The QD is defined by a harmonic external confin- 
ing potential v cxt (r) — uj 2 r 2 /2 on the xy plane, where lj is 
the confinement strength. The average electron density 
in the QD can be approximated by a density parameter 
r s = A^ 1 / 6 ^ -2 / 3 (Ref. [241). The parameter corresponds 
to the average radius of an electron in a QD with an av- 
erage number density uq = l/(7rr 2 ). In the case of the 
QS, the confining potential is a 2D rectangular quantum 
well with steep (hard-wall) boundaries ; 25 ' 26 and the den- 
sity parameter can be determined from r s = y/ A/ (irN), 
where A is the area of the QS. 

For both test systems, in the parameter ranges con- 
sidered here for N, u, and A, the LDA is known to 
yield very good total energies in comparison with exact or 
semi-exact many-electron methods, e.g., quantum Monte 
Carlo calculations. 25 ' 27 ' 28 Therefore we use the LDA en- 
ergies as the reference data in this work. We compute 
the LDA energies in the standard KS scheme by apply- 
ing octopus real-space DFT code^ For the LDA corre- 
lation [last term in Eq. (|26[) ] we use the parametrization 
of Attaccalite et al^ 

We apply our functional in Eq. (f2"3")l , as well as the TF 
expression in Eq. (|24|) . by using the self-consistent LDA 
density as the input density in a one-shot calculation. In 
this way, all the functionals are evaluated with the same 
particle densities, and the LDA densities may be consid- 
ered as reasonable estimations of the exact ones. The 
possibility for a self-consistent application of our func- 
tional, and its practical relevance, are subjects of future 
investigation. 
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FIG. 1: (color online) Relative error in the total energy 
calculated with our functional (filled circles) and with the 
Thomas-Fermi approximation (open circles) for a six-electron 
parabolic quantum dot as a function of the density parameter 



FIG. 2: (color online) Relative error in the total energy cal- 
culated with our functional (filled symbols) and with the 
Thomas-Fermi approximation (open symbols) for rectangu- 
lar quantum slabs with N = 12 . . . 200 as a function of the 
density parameter r s . 



IV. RESULTS 

First we consider the total energies of a spin- 
unpolarized quantum dot with a fixed number of elec- 
trons, N = 6, as a function of the density parameter 
r s . Figure Q] shows the relative error of our functional 
(filled circles) and the TF approximation (open circles) 
in the total energy against the reference LDA result, i.e., 
I -Elda Elba ■ We emphasize that for this system the 
relative error of the reference LDA is below 0.003 with re- 
spect to quantum Monte Carlo calculations.— The values 
used for r s correspond to a wide range of the interaction 
strength, covering well the typical values (r s ~ 1 ... 5) 
used when modeling QD's within the effective-mass ap- 
proximation at the interface of GaAs and AlGaAs.— 

Overall, we find an excellent agreement in the total en- 
ergies between the LDA and our functional. The relative 
error remains below ~ 5 % through the full range of r s . 
On the other hand, the TF approximation is accurate 
only close to the noninteracting limit (r s — > 0), whereas 
in general the TF error is dozens of percent. The over- 
estimation of the total energy in the TF approximation 
is plausible due to the lack of exchange and correlation 
energies which are both always negative. 

Next we focus on the rectangular QS and vary both r s 
and N. The results are given in Fig. [2] which, similarly 
to Fig. [TJ shows the relative total-energy errors of our 
functional (filled symbols) and the TF (open symbols). 
The electron number varies in the range N = 12 . . . 200. 
As in the case of a parabolic QD, the accuracy of our 
functionals is superior to that of the TF, except at small 
r s - 

A significant feature in Fig. [2] is the consistency of the 
accuracy of the present functional with N. Instead, the 
validity of the TF approximation strongly depends on 



the electron number, which is due to the fact that at 
fixed r s , the relative amount of (classical) Hartree en- 
ergy of the total energy increases with N. However, even 
at large electron numbers our functional is considerably 
more accurate than the TF approximation, on condition 
that the system is not too close to the noninteracting 
(small-r s ) regime. For example, at N = 120 and r s ~ 3.3, 
the relative errors of the present functional and the TF 
approximation are 2% and 10%, respectively. Hence, 
our functional is expected to be a reliable tool for total- 
energy calculations in systems that are computationally 
not reachable by, e.g., the LDA (see, e.g., Ref. Il5|). In 
fact, the N = 200 case in Fig. [2] was close to the numer- 
ical limit of our LDA calculations. 

We observed numerically that in QD's the difference 
between Ttf and Tks is very small. As it is known, it 
actually goes to zero in the limit r s — > (Ref. [l9h . In 
QS's, on the other hand, Ttf may largely underestimate 
Tks- However, the relative contribution of this under- 
estimation to the total energy strongly decreases as a 
function of N. It remains to be seen whether a fully self- 
consistent application of the presented functional may be 
carried out providing either accurate or at least better 
densities than the Thomas-Fermi ones. 



V. CONCLUSIONS 

In summary, we have derived an explicit density func- 
tional for the total energy of electrons in two dimensions. 
The functional is numerically highly efficient due to the 
freedom from orbitals and from the calculation of the 
Hartree integral. When applied to models of semicon- 
ductor quantum dots and slabs up to hundreds of elec- 
trons, and up to strong electron-electron interactions, we 
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have found a good overall agreement with respect to the 
local-density approximation, and a significant improve- 
ment over the Thomas-Fermi approximation. Natural 
future developments of the present work include the spin- 
dependent generalization and the capability to deal with 
dimensional crossovers (such as from two to three dimen- 
sions, or from two to one dimension). 
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